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We use semiclassical Hamiltonian optics to investigate the propagation of light rays through two- 
dimensional photonic crystals when slow spatial modulation of the lattice parameters induces mixed 
stable-chaotic ray dynamics. This modulation changes both the shape and frequency range of the 
allowed frequency bands, thereby bending the resulting semiclassical ray trajectories and confining 
them within particular regions of the crystal. The curved boundaries of these regions, combined with 
the bending of the orbits themselves, creates a hierarchy of stable and unstable chaotic trajectories 
in phase space. For certain lattice parameters and electromagnetic wave frequencies, islands of 
stable orbits act as a dynamical barrier, which separates the chaotic trajectories into two distinct 
regions of the crystal, thereby preventing the rays propagating through the structure. We show that 
changing the frequency of the electromagnetic wave strongly affects the distribution of stable and 
unstable orbits in both real and phase space. This switches the dynamical barriers on and off and 
thus modulates the transmission of rays through the crystal. We propose microwave analogues of 
the photonic crystals as a route to the experimental study of the transport effects that we predict. 

PACS numbers: 42.65.-k, 05.45.-a 



I. INTRODUCTION 

Quantum chaos explores the quantum-mechanical 
properties of systems whose classical counterparts ex- 
hibit deterministic chaos [TJ [2] . Most studies of quantum 
chaos in experimentally-accessible systems have focused 
on atoms [TJ [3J 0] and semiconductor systems, in which 
electron confinement by potential energy barriers of var- 
ious shapes generates chaotic classical trajectories. 

Semiconductor superlattices, comprising a series of 
quantum wells, have been used to realize a fundamentally 
different regime of nonlinear dynamics, in which the effec- 
tive classical Hamiltonian originates from an intrinsically 
quantum-mechanical feature of the system: tunnel cou- 
pling of the wells to form energy bands, known as mini- 
bands [5Hn|. When a bias voltage and tilted magnetic 
field are applied, the energy versus crystal momentum 
dispersion relation for the energy bands generates chaotic 
dynamics equivalent to those of a one-dimensional har- 
monic oscillator driven by a plane wave [5J [7HS] • Unusu- 
ally, this system does not obey the KAM theorem [ToT - 
I14j , which means that chaos switches on and off abruptly 
at critical values of the applied fields. At the critical 
fields, the electron phase space is threaded by intricate 
web patterns, known as stochastic webs [TUMl4| . These 
webs form a network of conduction channels, through 
which the electrons propagate in real space [5H9] . When 
the web is switched on, the electrons undergo chaotic 
diffusive motion along its filaments, thereby producing a 
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sharp increase in the measured and calculated direct cur- 
rent [71 [5] and increasing both the power and frequency 
of GHz/THz current oscillations by an order of magni- 
tude [5]. Consequently, chaotic semiclassical transport 
through energy bands provides a mechanism for control- 
ling transport through periodic structures. In the ab- 
sence of electron scattering, switching due to the for- 
mation and destruction of stochastic webs is extremely 
sensitive, due to the inherent instability of chaotic or- 
bits, and would produce S— function peaks in the current- 
voltage characteristics. In real superlattices, though, the 
peaks are broadened by electron scattering. This limita- 
tion can, in principle, be overcome by using an analogous 
system comprising ultracold atoms moving through the 
energy bands of an optical lattice with a tilted harmonic 
trap [15]. 

Due to the formal similarity between Maxwell's equa- 
tions and the Schrodinger equation pQ, microwave and 
optical systems have been used to study the ray-wave 
correspondence, analogous to quantum chaos, for elec- 
tromagnetic radiation. Much of this work has focused on 
electromagnetic billiards in which the radiation is con- 
fined in a quasi two-dimensional (2D) region of space, 
usually a cavity or dielectric block, whose boundary is 
shaped to create chaotic dynamics for light rays confined 
within the block. There has been considerable techno- 
logical interest in such systems because, when used as a 
resonator for a micro laser, they can produce strongly fo- 
cused emission whose power is orders of magnitude higher 
than obtained from traditional circular resonators jTBJ- 
[20] . Microwave systems have also been used to study the 
ray-wave correspondence in billiards, photonic crystals, 
and in a range of disordered dielectric media |2"1T[2"41 . 
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FIG. 1: Schematic diagram of the photonic crystal, which 
comprises two arrays of parallel dielectric sheets (gray) inter- 
secting at 90° . Each sheet is of infinite area. The sheets are 
parallel to the x-z or y-z plane. The x' and y' axes are related 
to the x and y axes by a rotation of 9 about the z axis (inset). 



In 2D photonic crystals [25, 26 , ray dynamics related 
to the chaotic semiclassical motion of electrons in a su- 
perlattice miniband has been predicted when the lattice 
parameters, and hence the local refractive index, vary 
slowly with position [27] [28]. Such structures offer major 
advantages over superlattices for studying, and poten- 
tially exploiting, chaotic band transport effects because 
they operate at room temperature and without an ap- 
plied magnetic field. As a ray passes through the crys- 
tal, changes in the lattice parameters alter the local pho- 
tonic band structure, in particular the shape and fre- 
quency range of the allowed bands. This exerts an ef- 
fective "force" on the ray [27] [29], whose trajectory can 
be determined by solving Hamilton's equations [57], just 
like calculating the path of a classical particle subject to 
real forces. When the lattice parameters are chosen to 
approximate the potential landscape experienced by an 
electron in a superlattice with a tilted magnetic field, the 
ray paths share many characteristics of non-KAM chaos. 
In particular, they have a rich mixed stable-chaotic phase 
space structure [27j . 

In this paper, we study the origin of mixed-stable- 
chaotic ray dynamics in spatially-modulated photonic 
crystals and how such dynamics depend on the param- 
eters both of the crystal and the electromagnetic wave. 
We focus on two key aspects of the ray dynamics. Firstly, 
we consider how the spatial variation of the lattice pa- 
rameters affects the range and dispersion of the allowed 
frequency bands and, hence, the ray trajectories obtained 
from Hamiltonian ray dynamics. Secondly, we analyze 
how the form of the ray trajectories, in real and phase 
space, depends on the frequency, u>, of the electromag- 
netic wave. In contrast to previous work |27j . which con- 
sidered ray dynamics at fixed frequency, we find that 
the location and extent of the stable and chaotic or- 
bits depend strongly on uj. By considering four critical 
points in the local Brillouin zone, we determine the loci of 
classically-allowed regions of the photonic crystal, where 
rays can propagate, and the corresponding forbidden re- 
gions that the rays cannot enter. Within the classically 
allowed regions, we identify additional dynamical barri- 



ers, which are formed by localized stable orbits and are 
impenetrable to chaotic ray trajectories [27]. We show 
that as the ray propagates through the crystal, changes 
in the local bandstructure cause it to reach the band 
extrema (highest or lowest frequency within an allowed 
band), where Bragg reflection occurs. Together, Bragg 
reflection and the action of the dynamical barriers limit 
the spatial region in which a ray moves, and drive many 
ray paths chaotic. As uj increases, a single region of chaos 
splits into two distinct chaotic seas separated by a dy- 
namical barrier. Further increasing u) suppresses chaos 
in one of these two regions, thereby deactivating the dy- 
namical barrier and its effect on ray transport. 

The structure of the paper is as follows. In Sec- 
tion [Tl] we define the crystal structure and consider the 
frequency-versus Bloch wavevector dispersion relations 
for the allowed frequency bands in which electromag- 
netic radiation can propagate through the crystal. Next, 
in Section III we introduce the semiclassical (Hamilto- 



nian) equations of ray motion and consider, qualitatively, 
how the shape of the dispersion relations influences the 
ray trajectories, in particular whether they are open or 



closed. In Section IV we consider the spatial regions of 
the crystal in which the rays propagate and show how the 
shapes of these regions, and the orbits themselves, are de- 
termined by the form of the dispersion relations and the 
photonic crystal structure. In Section[V] we demonstrate 
the existence and influence of dynamical barriers, which 
originate from islands of stability within a chaotic sea 
rather than from a classically forbidden region of high 
potential energy. In particular, we show that these dy- 
namical barrriers can control the transmission of electro- 
magnetic waves through the photonic crystal. In Section 
VI we discuss possible routes to the experimental real- 
ization and study of the dynamical regimes that we intro- 
duce. Finally, in Section |VII[ we summarize our results 
and draw conclusions. 



II. THE STRUCTURE OF THE CRYSTAL 

We consider photonic crystals created from two inter- 
secting sets of dielectric sheets (see Fig. [lj. One set is 
parallel to the x-z plane and the other is parallel to the 
y-z plane. The crystal is invariant, and taken to be infi- 
nite, along the z axis. This creates a photonic crystal in 
which the dielectric sheets enclose rectangular air gaps, 
as shown in Fig. [T] 

As in Refs. [35] [2jj], we characterize each dielectric 
sheet by the dielectric parameter m = (e r — l)d, where 
d and e r are, respectively, the width and relative per- 
mittivity of the sheet. At the lines where two sheets 
intersect, the relative permittivity is 2e r — 1. At position 
r = (x,y) within the crystal, we take the local unit cell 
length in the x and y directions to be l x and l y respec- 
tively. In the limit d— >0, if the electric field vector of an 
electromagnetic wave of frequency uj is parallel to the z 
axis, the local frequency versus wavevector dispersion re- 
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lation, ui\ oc (fj.,l x ,ly), of the wave satisfies the equations 

EH, 



cos (n x l x ) = cos (k x l x ) 



COS (fJ-yly) — COS (kyly) 



sin(A; x Z x ) (fa) 



'"■^loc. 

2c^ A^7/ 



sin (kyly), (fb) 



where c is the speed of light in vacuo, fi =(ii x ,/j, y ) is the 
Bloch wavevector, and k=(k x ,k y ) is the local wavevec- 
tor of the electromagnetic wave between the dielectric 
sheets. At time t, the local electric field magnitude, 
E = Eoe 1 ( k - r -"*) ) is determined by the wavevector k. 
The Bloch wavevector, /x, determines the spatial devel- 
opment of the phase of the electromagnetic wave as it 
propagates through the crystal [25] . 



When the y-component of the wavevector between the 
dielectric sheets is evanescent, so that k y = iq y where q y 



is real, Eq. ( lb ) can be written in the form (25j : 



cos (Hyly) = cosh (q y l y ) - 2 '° c sinh (q y l y ) . (2) 

Zc q y 




FIG. 2: Dispersion surface, CJioc(fJ>x, fly) calculated for l x = 
0.75 fim, l y — 0.5 pm. and mo = 2 /im. Three constant 
frequency contours are shown with lo = 5.7 x 10 14 rad s _1 
(dashed curve), lo = 6.75 x 10 14 rad s" 1 (solid curve) and 
lo — 7.35 x 10 14 rad s _1 (dot-dashed curve). At the bottom of 
the figure the contours are shown projected onto the cji oc = 
plane. Rays propagate perpendicular to these contours. One 
quadrant of the reduced Brillouin zone is shaded gray with its 
corners labelled A, B, C, D. The inset, top right, shows the 
repetition of the Brillouin zone and illustrates that the low- 
est and highest frequency contours (dashed and dot-dashed 
curves respectively) form closed curves while the intermedi- 
ate contour (solid) forms an open curve. 



If the lattice parameters m, l x and l y vary sufficiently 
slowly, the local dispersion relation at any point in the 
crystal is, to good approximation, the same as the disper- 
sion relation for an infinite crystal with a constant unit 
cell identical to the local unit cell [27]. Here, we choose 
m and l x values that vary slowly with position and are 
defined by the continuous functions l x — 1$ exp(— rjx) and 
m = TOq(1 — py' ), where the constants Z , m , r), and p 
are specified below. The x' — y' axes, are related to the 
x — y axes by a rotation of 9 about the z axis, as shown 
in Fig. 



The solutions of Eqs.(fal and ^ define a number of 
allowed frequency bands of the crystal. In the bands con- 
sidered here, we show below that changes in the positions 
of both the top and bottom of the band, i.e. the high- 
est and lowest frequencies, play a key role in driving the 
ray paths chaotic. For that reason, the lowest frequency 
band is of no interest here, because it always has a solu- 
tion extending down to lo = 0, which means that chaotic 
ray paths within the band are unbounded on the right 
(i.e. for high x) [27 ]. Instead, the band of interest in this 
work is the next lowest frequency band in which the rays 
have a locally evanescent component in the y direction 
and a locally propagating character in the x direction 

[IB HZ]. 



III. THE RAY TRACING 

The ray paths in the crystal are determined by the pair 
of Hamilton's equations [27J [3D] 



dt 

dfj, 

lit ' 



dH 
dr 



The 2D Hamiltonian is defined as [27] 

H(r, n, lo) = w loc (/x, l x , l y ) - lo, 



(3a) 



(3b) 



(4) 



where lo is the angular frequency of the electromagnetic 
wave and the local lattice constants l x ,l y both depend 
on position r. To determine ray paths in the second 
band we first specify lo and the initial position of the ray 
front, both spatially in the crystal, (x, y), and the Bloch 
wavevector in the band, (fi x ,p y ). We then solve Eqs. 



(3a) and (3b) numerically using a fourth-order Runge- 
Kutta method. It follows from Eqs. (3a I and Q that 



dr 

dt 



dLOloc 



(5) 
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From Eq.([5]) we see that the direction in which a ray 
propagates will always be perpendicular to the constant 
frequency contour u) = lu\ oc on the dispersion surface. 

Figure [2] shows ui\ oc versus fi x and fj, y at the point (x, y) 
— (28.7 /im, y — 11 /on), in the crystal where l x = 0.75 
/im, l y — 0.5 /im, and mo = 2 /im. We also take 9 = 21° 
and 77 = 10 m : parameters that will be used through- 
out this paper. The three constant frequency contours 
marked in Fig. [2] correspond to ui = 5.7 x 10 14 rad s _1 
(dashed curve), w = 6.75 x 10 14 rad s _1 (solid curve) and 
cj = 7.35 x 10 14 rad s _1 (dot-dashed curve). 

Due to the symmetry of the frequency band, we need 
only consider one quadrant of the reduced Brillouin zone 
(gray shaded region in Fig. [2]) for which < fi x < tt/I x 
and < [i y < ■k /ly. In fact, it is the four corners (marked 
A, B, C and D in Fig. [2]) of this reduced Brillouin zone 
quadrant i.e. the points where fi x = or ir/l x and fi y = 
or Tr/l y that will be of most interest here. 

The inset in Fig. [2]reveals that, since the Brillouin zone 
repeats with a period of 2n/l x in the \x x direction and 
lixjly in the /j, y direction, the lowest frequency contour 
at to = 6.75x 10 14 rad s _1 (dashed curve), and the highest 
frequency contour at ui — 7.35 x 10 14 rad s -1 (dot-dashed 
curve), both form closed curves. However, this is not the 
case for the intermediate contour at u> — 6.75 x 10 14 rad 
s _1 (solid curve), which is open in the extended zone 
scheme. 

Since the frequency band varies monotonically in both 
the \i x and /x y directions over the shaded section of the 
reduced Brillouin zone in Fig. [2j the gradient du>i oc /dfi is 
only parallel to the p, x axis along the lines fi y = or ir/l y , 
and only parallel with the /i y axis along the lines \i x = 
or tt/I x . When contours form open curves they are re- 
stricted to lie either between /i a = and ir/l yi or between 
fi, x = and tt/I x . Since a ray propagates perpendicular 
to the constant frequency contour where wi oc = to. if this 
contour forms an open curve there are directions in which 
the ray cannot propagate. For example, it can be seen 
from Fig. [2] that for the (solid) constant frequency con- 
tour at w = 6.75 x 10 14 rad s _1 , rays cannot propagate 
parallel, or almost parallel, to the y axis, because io\ oc is 
never parallel to the fi x axis. Below, we show that this 
restriction causes spatial localization of the rays in the y 
direction. 

For a given uj, a constant frequency contour forms a 
closed loop if u> = iO\ oc is between the bottom of the 
band, point A in Fig. [2] and the saddle point B. The 
contour is also closed if wi oc is between the top of the 
band, point D in Fig. [2] and the saddle point C. The 
constant frequency contours between the two frequencies 
corresponding to u>\ oc at points B and C always forms 
open curves. 



IV. THE SHAPE OF THE RAY PATHS 

Henceforth, we consider crystal parameters, L = 1 /im, 
l a = 1 /im, 8 = 21°, j] = 1 x 10 4 m _1 , m = 2 /im, 
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FIG. 3: Loci of points in the crystal where to — 6 x 10 1 rad 
s" 1 coincides with the uj\ oc values at points A (solid curve), 
B (dotted curve), C (dot-dash curve) and D (dashed curve) 
in the reduced Brillouin zone shown in Fig. [2] In the gray 
regions, u lies in a bandgap, and hence there are no propagat- 
ing solutions. Rays can only propagate in the white regions 
bounded by the solid and dashed curves where uj coincides 
with the band minimum and maximum, respectively, shown 
in Fig. [2] Within the black regions, m is negative, mean- 
ing that either e r or d must also be negative. Rays cannot 
reach these regions of negative permittivity because they are 
unable to traverse the outer gray area where to lies within a 
bandgap. Consequently, the regions of negative permittivity 
do not affect the ray dynamics that we consider. 



p = 1 X 10 10 m and, unless otherwise stated, take uj — 
6 x 10 14 rad s . When there is an evanescent component 
to the electric field, this is in the y direction. 

Spatial variation of the unit cell changes both the shape 
of the dispersion curve, LU\ oc (fi,l xl ly), and the range of 
frequencies that it spans. Thus, in parts of the crystal 
the electromagnetic wave may lie within an allowed band, 
but in other parts its frequency coincides with a bandgap. 

It is therefore essential to identify where in the crystal 
lu coincides with an extremum (points A and D in Fig. 
[2j or a saddle point (points B and C in Fig. [2]) in the 
dispersion surface. 

To illustrate this, Fig. [3] shows the four loci of points 
in the crystal at which w coincides with the two extrema 
or the two saddle points in the local u;i oc (/x, l x , l y ) surface 
(i.e. points A-D in Fig. |2|. The solid, dotted, dot-dashed 
and dashed curves correspond to points A, B, C and D 
respectively. Between curves A and D, the white region 
in Fig. [3j ui lies within the second allowed frequency 
band. Whilst it is not universally true that point B lies 
at a lower frequency than point C, as can be seen from 
the crossing of curves B and C at x = x cross ~ —48 /im in 
Fig. [3j this is the case throughout the region of interest 
here. The gray region in Fig. [3] shows the regions of the 
crystal for which there are no propagating solutions at 
frequency lu. In the black region, m < 0, meaning that 
either e r or d would need to be negative. Although some 




x (|um) 

FIG. 4: Ray paths calculated in a region of the crystal where 
all trajectories are stable. The four curves are the loci of 
points in the crystal where u — 6 x 10 14 rad s _1 coincides with 
the uui oc values at points A (solid curve), B (dotted curve), C 
(dot-dash curve) and D (dashed curve) in the reduced Bril- 
louin zone. These four curves strongly influence the shape of 
the ray paths and determine the turning points. 



materials do have negative values of permittivity, we do 
not consider them here. 

As discussed above, at any point in the crystal a ray 
can, in principle, propagate in any spatial direction pro- 
vided that the cl>i oc = uj contour forms a closed loop. By 
contrast, the range of possible propagation directions is 
limited for open contours. In Fig. [3j when x > x cross , at 
any point in the crystal that lies between curves A and 
B, or between curves C and D, the iso-frequency contour 
for u!\ oc = uj = 6 x 10 14 rad s^ 1 forms a closed loop. 
Conversely, for any point in the crystal between curves 
B and C, the iso-frequency contour forms an open curve, 
which restricts the range of directions in which a ray may 
propagate. 

Figure [4] shows a selection of stable ray paths in the 
crystal. The turning points of these rays are produced 
either by Bragg reflection, whenever \i x = ±ir/l x or fi y 
— ±7r/l y , or by the rays reaching band extrema away 
from the Brillouin zone boundaries as a result of the 
band changing shape throughout the crystal. Figure [4] 
reveals that the ray paths are shaped by the shifting of 
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FIG. 5: (a) Two chaotic ray paths propagating in two distinct 
regions of the photonic crystal when uj = 6.09 x 10 14 rad s _1 . 
The crystal parameters are l y = 1 lira, lo = 1 fim, 6 = 21°, 
r\ = 1 x 10 4 m _1 , m = 2 /im, p = 1 X 10 1() m. Panels (b) and 
(c) are enlargements of the left and right-hand boxes in (a). 
Crosses in (b) and (c) mark the points where Bragg reflections 
occur. 



the bandstructure, with curves A, B, C and D passing 
through extrema in the envelopes of the ray trajectories. 
Between curves A and B, Bragg reflections occur when 
fj, x = ±it/l x . At frequencies above that at point B in 
Fig. [2] such reflections no longer occur because the iso- 
frequency contours in Fig. [2] only intersect with the edge 
of the Brillouin zone where fi y = ±u /l y (between points 
B and D). 



V. DYNAMICAL BARRIERS 

Previous work has shown that the phase space of rays 
in the photonic crystal contains mixed stable-chaotic re- 
gions and that, for certain w, two or more chaotic regions 
are separated by islands of stability [27 . A trajectory 
starting in one chaotic sea cannot cross the stable island 
and is therefore unable to enter the second chaotic sea. 
Consequently, the stable island is known as a dynami- 
cal barrier. Although there is sometimes diffusion across 
dynamical barriers, such as via cantori (see, for exam- 
ple, [31]), in the photonic crystals considered here, the 
dynamical barrier can be completely impenetrable, lo- 
calizing chaotic rays to two distinct regions of the phase 
space. 
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In this section, we explain the existence of these dy- 
namical barriers and explore their effect on ray transport. 
The dynamical barriers occur for a wide range of uj and, 
at the higher frequencies, separate the chaotic regions in 
real space as well as in phase space. Consequently, dy- 
namical barriers provide a mechanism for controlling the 
transmission of light through photonic crystals. 

Figure [5]ja) shows two distinct chaotic paths in the 
photonic crystal that are almost separated in real space. 
Figures [5jb) and (c) show enlargements of the regions 
within the boxes in (a), in which the rays are driven 
chaotic. The solid curve, (A) , in Fig. [5] marks the loci 
of points where the bottom of the frequency band (point 
A in Fig. [2) is at w loc = 6.09 x 10 14 rad s -1 , whilst the 
dashed curve, (D), marks the loci of points where the top 
of the band (point D in Fig. [2]) is at cji oc = 6.09 x 10 14 
rad s _1 . 

The crosses in Fig. [5] (b,c) mark the points along the 
real-space ray trajectory where the Bloch wavevector, ft, 
is at an edge of the reduced Brillouin zone shown in Fig. 
[2] At such points, the electromagnetic wave Bragg re- 
flects along either the x or y direction, corresponding to 
a turning point of the classical ray. Most of these turn- 
ing points occur away from the frequency band extrema 
labelled A, B and D in Fig. [2] and therefore also away 
from the edges of the classically- allowed (white) region in 
Fig. [5](b,c). But, occasionally, Bragg reflection occurs 
at, or very near to, one of these extrema. For example, 
the crosses situated exactly on the dashed curve in Fig. [5] 
(b) , which is the locus of points in real space along which 
the electromagnetic wave frequency co coincides with the 
top of the frequency band at Point D in Fig. [2j corre- 
spond to Bragg reflections along the y direction occurring 
at Point D. 

The dashed curve in Fig. [5] (a,b) is also the bound- 
ary between regions where waves can (white) and cannot 
(gray) propagate. Note that for y < — 2 /im, this sepa- 
rates the allowed region into two segments that occupy 
distinct ranges of x values. As noted above, when rays 
approach the dashed curve in Fig. [5] they Bragg reflect. 
Since this boundary is so strongly curved, it makes the 
ray path highly sensitive to the initial conditions, thereby 
driving paths that reflect from it chaotic. 

These ray trajectories have some similarities with the 
classical orbits of particles in a Sinai billiard comprising 
a rectangular outer wall, which encloses a circular scat- 
tering barrier. Particles within the billiard move between 
these two impenetrable walls. Orbits that interact with 
the circular barrier are chaotic because reflections from 
that barrier couple the motion parallel and perpendicu- 
lar to the outer wall. Similarly, in the photonic crystal 
considered here, chaotic ray orbits are generated by re- 
flections near the dashed boundary labelled D in Fig. [5] 

By contrast, the trajectory of the right-hand path in 
Fig. [5] (a) and enlarged in part (c), is driven chaotic 
by reflection from the concave (right-hand) end of the 
white (classically-allowed) region. In Bunimovich sta- 
dia, comprising straight parallel edges with semicircular 
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FIG. 6: A Poincare section in which a point is plotted each 
time a ray turns from traveling along the negative to the posi- 
tive x-direction. The points are generated by three ray paths: 
the two chaotic trajectories shown in Fig. [5] and a stable 
ray path that separates them. Points produced by the upper 
path in Fig. [5] are marked with crosses, while those produced 
by the lower path are marked with circles. Points generated 
by a stable path in the dynamical barrier are marked with 
triangles. 



walls at either end, chaos occurs due to reflections from 
the curved ends, which mix motion parallel and perpen- 
dicular to the side walls. In a similar way, ray paths 
confined near the concave right-hand end of the white 
allowed region in Fig. [5] are also driven chaotic. 

Rays that do not encounter either the concave end of 
curve A or the convex region around barrier D remain 
stable, and hence form a dynamical barrier that separates 
the two chaotic paths in Fig. [5] 

Figure [6] is a Poincare section in which a point is plot- 
ted each time a ray turns from traveling in the nega- 
tive x-direction to the positive x-direction. Points cor- 
responding to three different ray paths are shown. The 
crosses are generated by the upper of the two ray paths 
in Fig. [5] (a), whilst the circles correspond to the lower 
ray path. The broad spread of these two sets of points 
in the Poincare section confirms that each ray is chaotic. 
The triangular points in Fig. [6] are generated by a stable 
ray path (not shown) forming the dynamical barrier that 
separates the two chaotic orbits in Fig. [5^a). Since this 
ray is stable, it generates points on a smooth continuous 
curve in the Poincare section. 

To fully understand the origin of the dynamical barrier, 
it is useful to consider in more detail the bandstructure 
changes, throughout the crystal, produced by the varia- 
tion of l x in the x direction, and m with y' . 

Figure [7] shows dispersion surfaces, to\ oc {p x , fi y ), and 
iso- frequency contours at, (a-c), three different spatial 
points in the crystal. It illustrates how variations in l x 
and m both change the bandstructure. The effect of 
changing l x can be seen by comparing Figs. pfa) and 
(b). In (a) x = /im, where l x — Iq, and in (b) x = 50 
/im, where l x = 0.607 Zq. In both cases, y' — 0, and 
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FIG. 7: Dispersion surfaces uj\ oc (fj, x , fi y ) calculated at three different positions in the photonic crystal: (a) x = /mi, y' = 
/im where l x — l y = 1 /mi and m = 2 /jm; (b) a; = 50 /im, y' = /mi where l x — 0.607 2o /jm, l y = 1 /im and m = 2 /im; (c) 
x = /mi, y' = 7.07 /im where 2a = 2 V = 1 /im and m = 1 /im. The iso-frequency contours, shown both on the dispersion 
surface and projected onto the ui\ oc — plane beneath, are equally spaced by 3 x 10 13 rad s _1 . As l x decreases, the band shifts 
to higher 6Ji oc , as can be seen by comparing the dispersion surface in (b) with that in (a). In addition, over large sections of 
the dispersion surface in (b), the iso-frequency countours are almost parallel to the fi y axis, meaning that rays are, in general, 
restricted to propagate in directions close to the x axis. Comparing (c) to (a) reveals that as m decreases, the band shifts to 
higher u}\ oc and also broadens. 



m = 2 /im. Reducing l x , i.e. moving from (a) to (b), 
shifts the band to a higher range of frequencies. In addi- 
tion, reducing l x makes wi oc far less sensitive to changes 
in fiy, as revealed by comparing Fig. [7^b) with Fig. [7^a). 

This reduced variation of o;i oc with fj, y means that far 
more iso-frequency contours form open curves in Fig. 
[7][b) than in Fig. [7]ja), and that the curvature of these 
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FIG. 8: A selection of rays in the crystal with u; = (a) 6.15 x 
10 14 rad s~\ (b) 6.87 x 10 14 rad s" 1 . Also shown are the 
loci of points in the crystal where ui — oj\ oc at points A (solid 
curve), B (dotted curve), C (dot-dash curve) and D (dashed 
curve) in the reduced Brillouin zone shown in Fig. [2] In the 
gray regions there are no propagating solutions at the given 

UJ. 



contours decreases. The reduction of the curvature of the 
contours restricts the range of spatial directions in which 
a ray may propagate. For example, in Fig. [7]Jb), most 
rays propagate in directions almost parallel to the x axis. 

The effect of changing m on uj\ oc (ijl Xi fj, y ) may be seen 
by comparing Figs. [7][a) and (c), for which m — 2 /im 
and m = 1 /im respectively with fixed l x = 1 /im. As 
m decreases, the band shifts to higher values of ui\ oc and 
broadens so increasing the range of frequencies for which 
closed contours are found about points A and D. 

We now use these changes in the bandstructure to ex- 
plain more fully the ray dynamics. Figure 15] shows rays in 
the crystal calculated for uj = (a) 6.15 x 10 44 rad s _1 and 
(b) 6.87 x 10 14 rad s _1 . The solid, dotted, dot-dashed 
and dashed curves show, respectively, the loci where wi oc 
coincides with points A, B, C, and D in the Brillouin zone 
(Figs. | and |7j. 

The first point to note is that the (white) region in 
Fig. [8^b) in which the rays may propagate is located 
at higher x values than in Fig. [8][a). This is because 
the band spans higher frequencies when l x is smaller (see 
Fig. [T]), i.e. when x is larger, and so at higher u the rays 
propagate further along the x axis. 

Comparing Fig. |8ja) and (b) reveals that at higher uj 
[Fig. [8[b)] each ray is more confined in the ^-direction. 
This is a consequence of the rays propagating at angles 
that are closer to the x axis at higher uj; an the effect 
that can be understood by comparing Figs. [7]ja) and 
(b). Hence, the higher the value of uj, the greater the 
restriction on the propagation angles, and the narrower 
the ray path along the y direction. At higher uj, the 
regions of the crystal in which rays can travel in any 
direction, i.e. the regions between curves A and B, and 
curves C and D, are reduced. This corresponds to the 
reduced variation of wi oc with fj, y at higher uj, as shown 
in Fig. p». 

If uj is sufficiently low, for this structure when uj < 
5.91 x 10 14 rad s _1 , rays in the crystal encounter both of 
the regions (i.e the concave and convex barriers in Fig. 
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[5]) that generate chaos. At higher to, the spread of each 
ray orbit in the y direction is reduced, and a ray that 
encounters one of the regions where it is driven chaotic 
does not reach the other. When to > 6.12 x 10 14 rad s -1 , 
a region forms in the crystal where rays will encounter 
neither of the strongly curved barriers. These rays are 
stable and collectively form the dynamical barrier. 

As to rises above ~ 6.6 x 10 14 rad s _1 , the two dis- 
tinct regions of the crystal that induce chaotic ray paths 
shrink, and become increasingly localized near the con- 
cave and convex barriers. This increases the separation 
of the two sets of chaotic orbits in real space. If to is in- 
creased further, the right-hand chaotic region disappears 
completely, as shown in Fig. [sjb) , meaning that the dy- 
namical barrier also vanishes. An examination of rays 
near the right-hand concave barrier [within the right- 
hand box in Fig. [5]ja)] reveals that, due to the shifting 
bandstructure throughout the crystal and the effective 
force that this exerts on the rays, such rays are more com- 
plex than those in Bunimovich stadia. Consequently, the 
exact reason for the suppression of chaos in this region, 
which occurs at high to, is not yet fully understood. How- 
ever, as the ray paths become increasingly restricted to 
a small region around the center of the concave barrier, 
chaos disappears because the barrier acts on the orbits as 
though it is locally flat. By contrast, the chaotic orbits 
created by the left-hand convex barrier in Fig. [5] persist 
for even higher to. 



VI. PROPOSED EXPERIMENTAL 
REALIZATION OF THE MODULATED 
PHOTONIC CRYSTAL AND ASSOCIATED RAY 
DYNAMICS 

Microwave analogues offer major advantages for in- 
vestigating the propagation of electromagnetic waves 
through dielectric media and the relation between the 
waves and underlying ray dynamics 1 . Due to their 
large scale, photonic crystals for microwaves can be made 
very accurately with imperfections far smaller than the 
wavelength of the microwaves. Such imperfections there- 
fore have a negligible effect on the wave transmission pro- 
cesses. For example, modern computer-controlled milling 
machines can produce holes in a Teflon dielectric block 
with a precision of ±0.05 mm: orders of magnitude 
smaller than the wavelength of microwaves. In addition, 
the electric field profile of the electromagnetic modes can 
be directly measured in microwave systems pQ, thus en- 
abling study of, for example, the link between modal form 
and the transmission coefficient of the system. 

Recent experiments on the transmission of microwaves 
through ID photonic crystals comprising Teflon sheets 
have revealed that absorption of the microwaves by the 
Teflon has little effect on the shape of the measured 
transmission spectrum, which is in excellent quantita- 
tive agreement with the corresponding calculations |24j . 
Since absorption does not destroy coherent wave phe- 



nomena, in particular band formation, we expect that it 
will have a similarly small effect on the band transport 
processes discussed above. Consequently, modulated mi- 
crowave photonic crystals fabricated by milling cm-scale 
air holes into a Teflon sheet, or high-refractive-index (3.3) 
ceramic-filled plastic |32j , offer a promising starting point 
for the experimental verification and study of the dynam- 
ical effects discussed above, and their effect on electro- 
magnetic wave transmission. We note that photonic crys- 
tals with a gradient index have been successfully realized 
experimentally for metamaterials [33 -38], and such sys- 
tems might also be suitable for experimental studies of 
band transport phenomena in modulated photonic crys- 
tals. 

Here, though, we focus on possible microwave real- 
izations of the modulated photonic crystals. To ensure 
maximal correspondence with the system shown in Fig. 
[5j and analyzed in detail above, the microwave pho- 
tonic crystal would have rectangular air holes separated 
by thin dielectric walls. The spatial variation of the 
hole and wall widths follows from the requirement that 

lx = l>o ex P(~ 7 l x ) an d m = (e r — l)d — m (l — py' 2 ), 
as discussed in Section [H] Providing the walls are much 
narrower than the air holes, the local dispersion rela- 
tion of the frequency bands will be accurately described 
by Shepherd's analytical model [S3] (see Section |ll|. Wc 
note, however, that exact dispersion relations can be de- 
termined for any unit cell geometry by using standard 
techniques such as expansion over plane waves. 

As in previous experiments [I] [5T] , quasi 2D microwave 
photonic crystals could be used because the width of the 
system along z does not significantly affect either the fre- 
quency band structure or the resulting ray dynamics. In 
such structures, the dielectric block is « 1 cm thick along 
the z direction and enclosed by copper plates, which form 
a microwave cavity. To study transport through the 2D 
photonic crystal, continuous microwaves would be in- 
jected through a wire antenna, of diameter « 0.1 mm, 
which passes through the top plate at a fixed position, r. 
A vector network analyzer would be used to measure the 
spectrum T(to,r,r') for transmission to a receiving wire 
antenna at a variable position r' in the bottom plate for 
to in the range 30 x 10 12 < to < 250 x 10 12 rad s" 1 [UEE]. 

A key advantage of using a vector network analyzer is 
that for each to value, the electric field profile, E(x,y), 
can be mapped experimentally by using a moveable bot- 
tom plate to scan the receiving wire across the structure 
and analyzing how T(to,r,r') changes with r' pQ. Since 
both the modulus and sign of E(x, y) can be determined, 
continuous waves can study pulse propagation via a sim- 
ple Fourier transform of T(lu, r,r') to the time domain 
|21) . This Fourier technique, which is simpler than gen- 
erating microwave pulses PQ , could be used to investigate 
the effects of dynamical barriers on pulse propagation by 
measuring how T(to, r , r') varies with the positions of the 
two antennae, with to, and with 9. 

For 9 values that generate a dynamical barrier, which 
separates chaotic trajectories into distinct regions of 
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phase space (Figs. [5] and [6]), the dynamical barrier will 
reduce the range of ray trajectories that directly link 
the two antennae, thereby producing a measurable de- 
crease in transmission between them. As shown in Fig. 
[HJ the location, and even the existence, of the dynamical 
barrier varies rapidly with u>. Consequently, we expect 
that T(oj, r, r') will also demonstrate high sensitivity to 
changes in uj. Providing the parameters of the photonic 
crystal are accurately known, we expect good quantita- 
tive agreement between experiment and theory, which 
will help to elucidate the effect of the dynamical barriers 
on ray and wave transmission. 

The aim of the proposed microwave experiments is 
to investigate the feasibility of using dynamical barri- 
ers to control electromagnetic wave propagation through 
spatially-modulated photonic crystals and so provide 
proof-of-principle demonstration of this concept. We em- 
phasize that dynamical barriers occur in many systems 
with a mixed stable-chaotic classical phase. In the elec- 
tromagnetic system that we consider here, their existence 
requires spatial modulation of the bandstructure, which 
couples motion in the x and y directions, but does not de- 
pend on the precise details of the lattice structure used 
to produce that modulation: for example the shape of 
the unit cell. Consequently, the dynamics that we con- 
sider have the potential for scaling to optical wavelengths, 
where the air holes are usually circular rather than rect- 
angular. A rectangular unit cell speeds up the ray calcu- 
lations, making it easier to explore the large parameter 
space available, but is not essential to generate dynami- 
cal barriers, which should occur in a range of modulated 
photonic crystal structures. 

VII. SUMMARY 

We have shown that spatial modulation of the lat- 
tice parameters and, hence, the local bandstructure of 
a photonic crystal affects the stability, form, and loca- 
tion of Hamiltonian ray paths within it. The forma- 
tion of convex or concave barriers at the edges of the 
classically-allowed region drives orbits that interact with 



those barriers chaotic. As uj increases, changes in the dis- 
persion relation for the second allowed band make the ray 
paths increasingly localized in the y direction. In turn, 
this confines the chaotic ray trajectories to smaller ar- 
eas around the convex and concave barriers at the edges 
of the classically-allowed region. The location and ex- 
tent of the chaotic trajectories depend strongly on the 
value of uj. For certain uj, only a single region of chaos 
exists. However, as uj increases, this splits into two dis- 
tinct regions of chaotic ray paths, located near either the 
convex or the concave barrier. These chaotic regions are 
separated by a dynamical barrier comprising stable ray 
orbits. When uj is just large enough for the dynamical 
barrier to form, this barrier is leaky because, in some 
regions of the crystal, chaotic orbits that are isolated in 
real space overlap in phase space. But as uj increases, 
the chaotic regions separate fully in both real and phase 
space. In this regime, the chaotic orbits cannot penetrate 
the dynamical barrier and are thus unable to propagate 
through the crystal. As uj increases further, the chaotic 
rays become so localized that those near the concave bar- 
rier disappear completely, thus removing the dynamical 
barrier. 

We conclude that rays in slowly-modulated photonic 
crystals exhibit unique and rich dynamics, which are 
highly sensitive to both the spatial structure of the crys- 
tal and the value of uj. Dynamical barriers formed by 
the stable rays provide a mechanism for controlling light 
transmission in certain frequency ranges. The exponen- 
tial separation rate of two chaotic rays with slightly dif- 
ferent initial conditions also strongly affects the propa- 
gation of the electromagnetic radiation. This mechanism 
could be exploited in optical switches or routers and may 
also enhance angular dispersion due to the superprism 
effect in photonic crystal structures [3jJ], which has po- 
tential for applications in wavelength division multiplex- 
ing devices. We hope that our results will stimulate the 
further work required to determine how ray chaos af- 
fects the transmission of electromagnetic waves through 
modulated photonic crystals, and to test the validity of 
Hamiltonian optics in this complex regime of nonlinear 
dynamics. 
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